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ABSTRACT 


Recently, ISahni, Shafieloo fc Starobinskv (12014!) combined two independent measurements 
of II{z) from BAO data with the value of the Hubble constant Hq = H{z = 0), in order to test 
the cosmological constant hypothesis by means of an improved version of the Om diagnostic. 
Their result indicated a considerable tension between observations and predictions of the ACDM 
model. However, such strong conclusion was based only on three measurements of H{z). This 
motivated us to repeat similar work on a larger sample. By using a comprehensive data set of 
29 H{z), we find that discrepancy indeed exists. Even though the value of D,m,oh^ inferred from 
Omh^ diagnostic depends on the way one chooses to make a su mmary statistics (weighted mean 
or the median), the persisting discrepancy supports the claims of ISahni. Shafieloo fc Starobinsl^ 


( 2014!) that ACDM model may not be the best description of our Universe. 


Subject headings: cosmology: observations - dark energy - methods: statistical 


1. Introduction 


The disco very of accelerat i ng expansion of 


where w can either be con stant - the 

198Sh ('a 


p = wp, 

so called “quintessence” (jPeebles et al 


las a 


the Universe ( Riess et al.lll998t iPerlmutter et al 


value characteristic fo r the fixed point attractor) 


1999ll brought us a mystery which became one 


of the most important challenges for modern cos¬ 
mology and theoretical physics. This phenomenon 
has been confirmed since then in manifold ways 
using different probes like supernovae la, acous- 


or evo l ving in time ( Chevalier fc Polarskil 12001 
Lindeil I20n3^ since the scalar field could be ex- 


tic peaks in the C MBR ( Bernardis et al. 2000t 
ISoergel et al. 2003 ), bar^ Aon acousti c oscillations 
(BAO) ( Eisenstein et al.i2005 . 2011 ) in the distri¬ 
bution of the large scale structure etc. By now all 
this observational evidence was concordant with 
the simplest assumption that there exists non¬ 
vanishing cosmological constant A. While being 
the simplest, such a model is not theoretically 
satisfactory. There is a huge discrepancy if one 
tries to motivate A as a zero-point quantum vac¬ 
uum energy. Therefore alternative explanations 
appeared invoking the s calar field which settle d 
down in an attractor ( Ratra fc PeeblesI 1988[) . 
This motivated to push forward the phenomeno¬ 
logical picture of the so called dark energy de¬ 
scribed as a fluid with barotropic equation of state 


pected to evolve in time. The main drawback 
of such an approach is that it makes explicit as¬ 
sumption about the dark energy before its spe¬ 
cific model (a quintessence or evolving equation of 
state) can be tested on observational data. More¬ 
over, altern atives to dark energy such like modi¬ 
fied gravity (Dvali et al.ll2 0^ Sotiriou et al.ll20ir)t 


Noiiri et'Sl 20 lit Bengochea et al.l ~09[) cannot 

be easily tested within this phenomenology. All 
observational tests of quintessence pin-point its 
value close to ic = — 1 (within the uncertainties) 
which is equivalent to the cosmological constant. 
On the other hand tests of varying in time cosmic 
equation of state are much less restrictive and do 
not allow to make a decisive statement whether 
dark energy equation of state evolved or not. 

Therefore, we clearly need alternative probes: 
capable to discriminate between the cosmologi¬ 
cal constant and evolving dark energy not rely- 
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ing on the dark energy assumption and its equa¬ 
tion of state parametrization. One promising 
approach to such probes has been i nitiated by 
Sahni. Shafieloo fc Starob inskv ( 2008li and devel¬ 


oped further in (IShafieloo. Sahni fc Starobinskv I 
2012r) . By properly rearranging the equation for 
the Hubble function in the flat ACDM model: 
H'^{z) = -b -b 1 - Hm], they no¬ 

ticed that the so called Om{z) diagnostic (where 
h = H{z)/Hoy. 


Om{z) = 


h^{z) - 1 
(l + z)3-l 


( 1 ) 


should be constant and equal exactly to the 
present mass density parameter, if the ACDM 
model is the true one: Om{z)\cDM = This 

is remarkable and differentiates between ACDM 
and other dark energy models (including evolv¬ 
ing dark energy). Let us remark that essen- 
tially the same idea has also been formulated by 


Zunckel and Clarkson (|2008ll who called it “a lit¬ 


mus test” for the canonical ACDM model. Devel¬ 


oping this method (IShafieloo. Sahni fc Starobinskv 


oping 

1201211 they also considered a generalized two-point 
diagnostics 


Oto(zi,Z2) = 


ft^(zi) - h'^{z2) 

(1 + Zi)3-(1 + Z2)3 


( 2 ) 


which should also be equal to within the 

ACDM model but has an advantage that having 
H{z) measurements at n different redshifts, one 
has two point diagnostics, hence a consid¬ 

erably increased sample for inference. 

In t heir latest oaoer lSahni. Shafieloo fc Starobinskv 
( 2014ll used three accurately measured values of 
H{z) to perform this test. These were: the 


Hjz = 0) me asurement by (|Riess et al.l 12011 


lAd^ et al.ll2013), H(z = 0.57 ) meas urement from 
SDDS DR9 ( Samushia et al.l [201^ and the lat¬ 
est H{z) = 2.34 me asurement from the L ya for¬ 
est in SDSS DRll ( Delubac et al. 2015 1. They 
found that all three values of the two-point di¬ 
agnostics Om{zi,Z 2 )hy were in stro ng tension 
with the reported b y Planck (lAde et al 


could be of a paramount importance for the dark 
energy studies, an update of this test with a larger 
sample of H{z) is essential. 

2. Data, Methodology and Results 

As a basic dataset we used a sample of 29 
H{z) measurement s take n from the compila¬ 
tion of IChen et'H! ( 2014 1 modified in the fol¬ 
lowing wa y: one data point at z = 0.6 com¬ 


ing from Blake et al. (201^ 


was 


add e d and 

two data points from iCatzanaga et al. I ( 2009h 


were withdrawn. The reason for deleting afore¬ 
mentioned two points is that these results haye 
been debated in subs e quent papers e.g. by 


Miralda-Escude 


subse quent papers e.g. 
( 200911. kazin et al. I (l2010ll 


or 


Cabre and Gatzanaga (12 011 ). For the sake of 


consis tency with ISahni. Shafieloo &: Starobinsky 
(|2014ll . we haye also taken the la test BAO mea¬ 
surement by Delubac et al.l ( 2015h H{z = 2.34) = 
22 2 ± 7 instea d of H jz = 2.3) = 224 ± 8 
of Bnsca et al.l ( 2013l l. After these changes 
our data are essentiall y like the ones u s ed by 
Farooa &: Ratra I ( 2013ll with Busca et"^ ( 20Ly 
measurement replaced by IPelubac et al.l ( 2015ll . 
During the analysis we also made assessments 
on subsamples of this biggest one, as we will ex¬ 
plain further. Part of the data comes from cosmic 
chronometers - sp ectroscopy of galaxies a ssumed 
to eyolye passiyely ( Jimenez fc Loeb1l2002 ). Here¬ 
after, this differential age approach will be quoted 
as DA for short. The other part of the data comes 
from BAO - including the data poin ts used by 
Sahni. Shafieloo fc Starobinsky ( 2014h . Data are 
summarized in Tabl e 1. Then we proceed in ex¬ 
actly t he same way as ISahni. Shafieloo fc Starobinsky 
(I 2 OI 4 I I. i.e. for each pair of redshifts {zi,Zj) we 
calculate the improyed Om diagnostic: 


Omh^{zi, Zj) = 


h^jzi) - h^(zj) 

(1 + z,)3 - (1 + z,r 


(3) 


where: h{z) = H{z)llOQkmlseclMpc is di¬ 
mensionless Hubble parameter. In particular 
case of ACDM model the improyed diagnostic 
Eq. (|3]) should be equal to which luckily 

is the quantity best cons trained by the CMBR 


witn tne iimfin reportea b y riancx iiAde et al.l is tne quantity best cons trainea by tne 

2013 ). It has been noticed IPelubac et aT 20I5I: data, e.g. from Planck ( Ade et al. 20131) . Be 


Sahni. Shafieloo fc Starobinsky 20141) ttiat such 


result could be in tension not only with the ACDM 
model but with other dark energy models based 
on the General Relatiyity. Because this conclusion 


cause the sample of 29 data points leads to 406 
different pairs, we summarize our calculations 
on Figure [TJ One can see that the distribu- 

^m,0^ 


tion of inferred 0^ nh'^ is skewed and centered 
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around the value different from this r eported by 
Sahni. Shafieloo fc Starobinskv ( 2ni4l i. 


If one is to make a summary statistics one can 
do it in two ways. First, straightforward way 
would be to calculate the weighted mean: 


ESi Omh^z^, 


Omh 


{w.m.) 


1 ^/^Omh^ ,ij 


and the standard deviation: 


29 


- 1/2 


lEE'/ ^Omh^ ,ij 

,2=1 j>i 


(4) 


(5) 


where: 


^Omh^ ,ij 


((1 + Z .)^-(1 + Z ,) 3)2 


( 6 ) 


and (Jh^zi) denotes the uncertainty of the i — th 
measurement. We also assume that redshifts are 
measured accurately. This well known and often 
used approach relies, however, on several strong 
assumptions: statistical independence of the data, 
no systematic effects, Gaussian distribution of the 
errors. These assumptions, especially the Gaus- 
sianity of errors are not valid here. Hence the 
weighted mean of 


Omh 


{w.m.) 


= 0.1253 ±0.0021 


(7) 


is not a reliable measure, as one can see form the 
histogram on Figure [T] We will comment more on 
this non-Gaussianity in a moment. 

Therefore we took another, much more robust 
approach: to calculate the medi an. This approach 
was pioneered by the paper of Gott et id] ( 200lh 
and then used by the others , e.g. quite recently 
by (|Crandall &: Ratra 1120121) . The robustness of 
this method stems from the fact that if systematic 
effects are absent, half of the data is expected to 
be higher and another half - lower than the me¬ 
dian. Then, as the number of measurements N 
increases, calculated median approaches its true 
value. So the median has clear and robust meaning 
without need to assume anything about the error 
distribution. From the definition of the median, 
probability that any particular measurement, one 
of N independent measurements is higher than the 
true median is 50%. Consequently, the probability 


that n observations out of the total number of N 
is higher than the median follows the binomial dis¬ 
tribution: P = 2~^N\/[n\{N — ny]. This allows 
to calculate in a simple manner the confidence re¬ 
gions (e.g. 68% confidence region) of the median 
value estimated from the sample. Proceeding this 
way we have obtained: 


Ornh‘‘ _ n 1 c;c:n+0.0065 

{median) — QQ72 


( 8 ) 


In order to facilitate comparison between the 
inferred values of obtained from two statis¬ 

tical approaches and the Planck data, we display 
the results in Figure [5] 

Inconsistency between Omh^{zi, Z 2 ) diagnostic 
calculated on H{z) data and the Planck value of 
O,m,oh^ as well as mutual inconsistency between 
weighted averaging and median statistics schemes, 
motivated us to make some more detailed tests. 
First, we recalculated Omh^ for three sub - sam¬ 
ples: excluding the highest redshift z = 2.34 mea¬ 
surement, using only DA data and using only BAG 
data. Results are visualized on Figure[2]and shown 
in more details in Table 2. One can see that 
z = 2.34 point had a big leverage on the weighted 
average value — dropping this point one achieves 
agreement with the ACDM Planck value. How¬ 
ever, the question remains whether the weighted 
average scheme is approp riate. Therefore, fol¬ 
l owing Chen et al. ( 2003ll and ICrandall et al. I 
(|2014r) . we have drawn histograms of distribution 
of our measurements as a function of the num¬ 
ber of standard deviations away from central 
estimates (weighted mean and the median respec¬ 
tively) . Because of limited space we do not show 
them here, but report in Table 2 corresponding 
percentage of the distribution falling within ±1 ct 
i.e. liVo-l < 1. One clearly sees that they strongly 
deviate from the Gaussian 68 % expectation. We 
also tested the N„ distribution with Kolmogorov- 
Smirnov test which strongly rejected the hypoth¬ 
esis of Gaussianity in each sub-sample (with p- 
values ranging from 10“'^ to 10“^). Therefore 
our conclusion is that weighted average scheme is 
not appropriate here and the median statistics is 
more reliable. Both statistical methods: weighted 
mean and median produce similar results for the 
BAO data, but with addition of DA data these 
two schemes give drastically different results. This 
may suggest the existence of some systematic error 
in DA data. It is not obvious by itself, because the 
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Table 1 


Data of the Hubble parameter H{z) versus the redshift z, where H{z) and (jh are in km 
s“^ Mpc“^ These are essentially the data of IFarooq fc Ratra ( 2013), with the BA O 
measurement at the largest redshift H{z = 2.34) TAKEN AFTER [DELUBAC ET AL.I (j2015h . 


H{z) 


cth 


Method 


0.07 

69 

19.6 

DA 

0.1 

69 

12 

DA 

0.12 

68.6 

26.2 

DA 

0.17 

83 

8 

DA 

0.179 

75 

4 

DA 

0.199 

75 

5 

DA 

0.2 

72.9 

29.6 

DA 

0.27 

77 

14 

DA 

0.28 

88.8 

36.6 

DA 

0.35 

82.7 

8.4 

BAO 

0.352 

83 

14 

DA 

0.4 

95 

17 

DA 

0.44 

82.6 

7.8 

BAO 

0.48 

97 

62 

DA 

0.57 

92.9 

7.8 

BAO 

0.593 

104 

13 

DA 

0.6 

87.9 

6.1 

BAO 

0.68 

92 

8 

DA 

0.73 

97.3 

7 

BAO 

0.781 

105 

12 

DA 

0.875 

125 

17 

DA 

0.88 

90 

40 

DA 

0.9 

117 

23 

DA 

1.037 

154 

20 

DA 

1.3 

168 

17 

DA 

1.43 

177 

18 

DA 

1.53 

140 

14 

DA 

1.75 

202 

40 

DA 

2.34 

222 

7 

BAO 
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Fig. 1.— Histogram of Omh^ calculated from pairwise Om{zi, Z 2 )h'^ statistics performed on a sample of 29 
measurements of H{z) based on differential ages of galaxies and BAO. 


nonlinear relation between input variables {H{z)) 
underlying the Omh^ diagnostic may be the source 
of asymmetric uncertainties of the latter. How¬ 
ever, the fact that BAO and DA data deviate from 
the ACDM expected result in opposite directions 
strongly supports suggestion of unaccounted sys- 
tematics. This will be the subject of a separate 
study. 


3. Conclusions 


In this paper we attempted to assess the 
Omh?{z^ ,Z 2 ) diagnostic introduced and deve l- 


oped by IShafieloo. Sahni fc Starobinskv (2012 h . 


Th e main reason to do so was the recen t paper 
by ISahni. Shafieloo &: Starobinskv ( 2014 ) where 
they claimed that recent precise measurements of 
expansion rates at different redshifts suggest a se¬ 
vere tension with the ACDM model. We repeated 
this on a much more comprehensive data set of 29 
H{z) obtained by two techniques: DA and BAO. 
One can see from the Table 1 that uncertainties 
of H{z) obtained by different methods are dif¬ 
ferent. Even within the same methodology (DA) 
uncertainties are different from case to case. The 


Omh?'{zi, Z 2 ) diagnostic, involving ratio of certain 
differences (see Eq. ([3|)) calculated on our data 
produces an asymmetric distribution. This means 
that the weighted mean is not a reliable summary 
measure. Therefore we used a more robust ap¬ 
proach to calculate the median. 

Our result is that the value of Omh^ inferred 
from Omh^ diagnostic is indeed in tension with 
the one obtained by Planck (under assumption of 
the ACDM mode l). In our case, this tension is not 


so sev ere as m iSahni. Shafieloo fc Starobinskv 
(l2014ll (Om/i2 = 0.122 ±0.01 vs. ^mAhanck = 
0.1426 ± 0.0025). Even though the inferred value 
is sensitive to the way one chooses to make a sum¬ 
mary statistics: the weighted mean value is lower 
and the median value is higher than that obtained 
by Planck, they are both discrepant with each 
other. Non-Gaussianity in the data suggests that 
median statistics approach is more appropriate, so 
this tension cannot be alleviated by excluding high 
redshift data. 


Th is supports the claims of lSahni. Shafieloo fc Starobinskv 
(2 o 3) that the concordance model (ACDM) 
might not be the true or even the best one describ- 
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Fig. 2.— Values of Omh^ calculated from the Omh^ diagnostic as a weighted mean of data and as the median 
value. The result from Planck experiment is shown for comparison. Bands display the 68% confidence regions. 
Upper - left figure corresponds to the full sample n = 29 data points, upper - right shows the results when 
z = 2.34 point was dropped {n = 28), lower - left corresponds to DA data only (n = 23) and lower - right 
figure — BAO only (n = 6). 


Table 2 

Values of OmK^ diagnostic central values (weighted mean and median) and their 
“NON-GAUSSIANITY” indicated by the percentage of distribution falling within \N^\ < 1 FOR 
THE MAIN SAMPLE AND DIFFERENT SUB-SAMPLES. OBSERVATIONS OF THE CMB FROM PLACNK 

INFORM US THAT ^mfih^pianck = 0.1426 ± 0.0025. 


Sample 


\Ncr{w.m.)\ < 1 


\Na{median)\ < 1 

Full sample (n = 29) 

0.1253 ±0.0021 

80.54% 

0.1550lO:™f2 

75.62% 

z = 2.34 excluded (n = 28) 

0.1404 ±0.0047 

77.78% 

0-16821°™^^ 

82.80% 

DA only (n = 23) 

0.1448 ±0.0057 

77.47% 

0 1 Qc:0+0.0032 

U.iOOZ_Q QQ7g 

86.56% 

BAO only (n = 6) 

0.1231 ±0.0045 

100 % 

0 191 o+O.OOll 
U.iZi5 Q QQ35 

100 % 
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ing our Universe. Therefore, we also performed 
a quick test whether XCDM or CPL models 
(simplest evolving equation of state parametriza- 
tion) best fitted to the Planck or WMAP9 
data agree better with H{z) data. In such 
a case the Omh'^{zi,Zj) diagnostic defined be 
Eq.dll) will no longer be a single number 
but rather: Omh^{zi,Zj) — + (1 “ 

for XCDM and 

Omh'^{zi,Zj) = + (1 - 

— —SwgZj 

CPL parametrization. Therefore, for each pair 
we formed the residuals R{zi,Zj) by subtract¬ 
ing the right-hand sides from the left-hand sides. 
Because such residuals inherit non-Gaussianity 
from Omh^{zi, Zj) we summarized our findings 
with weighted average R(w.m.) and the median 
R(median) ■ If & particular model (XCDM or 
CPL) agreed better with the H{z) data than 
ACDM, then its R should have been closer to 
zero than R{kCDM)(^^ Tn '^ = —0.0173 ± 0.0033 
or R{KCDM)(^,d^au) = 0.0124t°;°o76- ^ the 
XCDM model we too k the w = --1.0507 ~*~n ntn? pa¬ 
rameter according to Cai et al. ( 2014ll best fit to 
Plack-|-WMAP9 data. For the CPL parametriza¬ 
tion of the equation of state, we used the val¬ 
ues: Wo = = 0.35to;49 best fitted 

to WMAP-|-CMB-|-BA O-H7n-|-SNe according to 
Hinshaw et al. I ()2013ll . 


The result is: R{XCDM)(^^ ^n.) = —0.0176 ± 
0.0025 and R{XCDM)^^,d,an) = 0.015ll°°“3^ 
For the CPL varying equation of state parametriza¬ 
tion we get R{CPL)i^y^ = —0.0275 ± 0.0073 or 
R{CPL)^^,d^an) = -0.05171° We see that 
they do not reconcile the discrepancy but their 
performance is even worse than ACDM. However, 
it is not a decisive conclusion, because what re¬ 
mains to be done is find the values of equation of 
state parameters w or {wo,Wa) best htted to the 
H{z) data (according to Omh^ diagnostics). This 
will be a subject of a separate study. 


Having confirmed the discrepancy between 
ACDM and H{z) data its origin should be stud¬ 
ied in greater details. One reason could be that 
our phenomenological description of accelerated 
expansion with ACDM is incorrect. On the other 
hand we pointed out that the conclusion (more 
specifically the direction of this discrepancy) de¬ 


pends on the statistical approach taken. So one 
should investigate possible systematics in both 
methods — DA and BAO and their effect on the 
conclusion. This is a subject of an ongoing study. 
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